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We present a gauge-invariant formalism to study the evolution of curvature perturbations in a 
Friedmann-Robertson-Walker universe filled by multiple interacting fluids. We resolve arbitrary 
perturbations into adiabatic and entropy components and derive their coupled evolution equations. 

We demonstrate that perturbations obeying a generalised adiabatic condition remain adiabatic in 
the large-scale limit, even when one includes energy transfer between fluids. As a specihc application 
we study the recently proposed curvaton model, in which the curvaton decays into radiation. We 
use the coupled evolution equations to show how an initial isocurvature perturbation in the curvaton 
gives rise to an adiabatic curvature perturbation after the curvaton decays. 
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I. INTRODUCTION 

The primordial curvature perturbation plays a central role in modern cosmology. It characterises large-scale density 
perturbations in our Universe from which smaller scale structures form via gravitational instability. Therefore much 
effort has been devoted to understanding the evolution of the curvature perturbation on large-scales in a general 
cosmology. A gauge-invariant formalism for cosmological metric perturbations was developed W Bardeen Q and the 
curvature perturbation (on uniform density hypersurfaces) C was introduced by Bardeen et al. [jildl shortly afterwards 
as a convenient gauge-invariant variable which remains constant for purely adiabatic perturbations on large scales. 
On large scales in an expanding universe it is essentially equivalent to the comoving density perturbation Hii 

The constancy of the curvature perturbation C in the case of a single perfect fluid follows directly from the local 
conservation of the energy-momentum tensor, in a suitably defined larce-scale limit 0- But C can change on arbitrarily 
large scales due to a non-adiabatic pressure perturbation 0, 0, IE 01 • Thus in a multi-fluid system it is in general 
necessary to follow the coupled evolution of curvature and entropy (or isocurvature) perturbations in order to determine 
the late-time curvature perturbation. 

There has been increasing interest in multi-field inflationary models and the spectrum of curvature ilHIIllIl 
IT^ and isocurvature [TE [lE| perturbations that may be produced and their correlations min In particular it 
has recently been suggested that the large-scale curvature perturbation C may be g enerated by initial isocurvature 
perturbations in a “curvaton” held which subsequently decays into radiation [l^ [l^ . 

Kodama and Sasaki developed a general formalism to describe the evolution of cosmological perturbations 
with multiple fluids (with corrections given in |2,‘ll). This formalism has subsequently been used by a number of 
authors 1^ 1^ 1^ |23 (see also 1^ I3(T| ). In particular Kodama and Sasaki applied their formalism to a 
matter-radiation huid in |24 1^ , where energy transfer can be neglected. One can argue on general grounds 0 that 
the entropy perturbations evolve independently of the curvature perturbation on large scales, but that the evolution 
of the large-scale curvature is sourced by entropy perturbations. Nonetheless there has been no detailed study of the 
evolution in general of curvature and entropy perturbations including energy transfer. By contrast a formalism to 
study the coupled evolution equations for curvature and entropy perturbations in models with multiple interac ting 
scalar helds has recently been developed by Gordon et al. and applied in a variety of scenarios 112,112,13 
(see also |12llEl)- 

In this paper we introduce a gauge-invariant formalism to follow the coupled evolution of curvature and entropy 
perturbations in multi-huid cosmologies when energy transfer between fluids is included. As an example we study 
the evolution of curvature and entropy perturbations in a curvaton scenario where the decay of the curvaton field 
represents the transfer of energy from the curvaton to radiation. We compare the results of numerical solutions of 
the coupled equations with analytic estimates based on the sudden decay approximation, where the curvaton and 
radiation are assumed to be non-interacting up until a given decay time. 
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II. GOVERNING EQUATIONS 

In this section we give the governing equations for the general case of an arbitrary number of interacting fluids in 
general relativity. We will consider linear perturbations about a spatially-flat FRW background model, as defined by 
the line element 


= —(1 + + 2aB^idtdx^ + a? [(1 — 2'tlj)Sij + 2E^ij] dx^dx^ , ( 1 ) 

where we use the notation of Ref.0 for the gauge-dependent curvature perturbation, ip, the lapse function, (p, and 
scalar shear, y = apE — aB. 

Each fluid has an energy-momentum tensor The total energy momentum tensor covariantly 

conserved, but we allow for energy transfer between the fluids, 

=Qu , (2) 

where is the energy-momentum transfer to the a-fluid, which is subject to the constraint 

EQ(«)= 0 - ( 3 ) 

Ol 

The equations hold for any type of fluid, the only requirement being the local conservation of the total energy- 
momentum tensor, = 0. 


A. Background equations 


The evolution of the background FRW universe is governed by the Friedmann constraint 

8ttG 


= 


-P, 


and the continuity equation 


p=-3H{p + P) , 


( 4 ) 


( 5 ) 


where the dot denotes a derivative with respect to coordinate time t, El = d/a is the Hubble parameter, and p and P 
are the total energy density and the total pressure 


'y Pa — P ^ 'y ^ Pa — P ■ 


( 6 ) 


The continuity equation for each individual fluid in the background is thus ^2 

Pa — 3-H (Pa T Pa^ “f Qa ; 


( 7 ) 


where the energy transfer to the a fluid is given by the time component of the energy-momentum transfer vector 
Q(a) ~ *3“ background. Equation implies that the background energy transfer obeys the constraint 


J2Qa=0. 


( 8 ) 


B. Perturbed equations 

Perturbing the constraint equation 0 yields the first-order equation [iilil 

(^ip + H(p^ — ii’ + Px) = —4:7rG6p , (9) 

where the comoving spatial Laplacian is denoted by jdx'^'^, and the momentum constraint equation (identically 

zero in the FRW background) is given by [a [23, [33 


ip + H<p = —AirGdq , 


( 10 ) 
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where <5/9 is the density perturbation and Sq the scalar S-momentum potential. 

Perturbing the continuity equation © yields an evolution equation for the total density perturbation 1^ 

6p + 3H{6p + 6P) - {p + P) i-tp - T + (P + P)x] = 0 ^ 


while total momentum conservation is given by MM 
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Sq + msq + {p + p)(i3 + 5P + n = 0, 

3 

where II is the total anisotropic stress. 

The perturbed energy transfer vector, Eq. including terms up to first order, is written as 

Q(a)0 S- SQ a 5 

Qc 




and Eq. © implies that the perturbed energy and momentum transfer obey the constraints 

^/a=0. 

Ct Ct 

The perturbed energy conservation equation for a particular fluid, including energy transfer, is then given by 

Sp^ + ^H{Spa + SPa) — {pa + Pa) ‘S'4’ H- T \^^ol + [Pa + ^ct)x] = Qa(p + SQa , 


while the momentum conservation equation is 


Sqa + iHSqa + {Pa + Pa)4’ + SP^ + 


2 


TIa — Qo 


Sq 


3^2— -"“p + p 


+ fa 


( 11 ) 

( 12 ) 

( 13 ) 

( 14 ) 

( 15 ) 

( 16 ) 

( 17 ) 


where the density, pressure, momentum and anisotropic stress perturbations of the individual fluids are related to the 
total density, pressure, momentum and anisotropic stress perturbations by 


'^Spa = Sp , SPa = SP , Y^qa = Sq , Ila = n . 


( 18 ) 


C. Gauge-invariant perturbations 

Both the density perturbation, Spa, and the curvature perturbation, ip, are in general gauge-dependent. Specifically 
they depend upon the chosen time-slicing in an inhomogeneous universe. However a gauge-invariant combination can 
be constructed which describes the density perturbation on uniform curvature slices or, equivalently the curvature of 
uniform density slices. 

The curvature perturbation on uniform total density hypersurfaces, Q is given by 


p 

while the curvature perturbation on uniform a-fluid density hypersurfaces, is defined asQ 

Ca = -V- . 

Pol 

The total curvature perturbation ca is thus a weighted sum of the individual perturbations 


( 19 ) 


( 20 ) 


( 21 ) 
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while the difference between any two curvature perturbations describes a relative entropy (or isocurvature) perturba¬ 
tion 

= 3(C„ - C/3) = -3i/ . (22) 

\ Pa P0 J 

The classic example of just such a relative entropy perturbation is a perturbation in the primordial baryon-photon 
ratio (with negligible energy transfer between the two fluids) 

Sb^ = 3(Cb - C 7 ) = - T— • (23) 

Pb 4 py 

This is also described as an initial isocurvature baryon density perturbation as Ss-y 5 pb/pb in the limit pb!P y —> 0. 

From the definitions of the total curvature perturbation m and the entropy perturbation we get 

Ca = C + o ^ ■ (24) 

a P 


D. Long-wavelength limit 

To describe the evolution of long-wavelength perturbations we will work in the ‘separate universes’ picture 0 where, 
smoothing over sufficiently large scales, the universe looks locally like an unperturbed (FRW) cosmology. Specifically 
we assume that we can neglect the divergence of the momenta in the zero-shear gauge, V‘^{5qa + {pa + Pa)x)i in 
Eq. lltill . 

In this long-wavelength limit, the perturbed continuity equation m becomes 

6p + m{5p -b (5P) = 3 (p -b P) . (25) 

Re-writing this equation in terms of the total curvature perturbation, C in Eq. m, gives Hi 

C =-, (26) 

p -b P 

where the non-adiabatic pressure perturbation is <5Pnad = SP — c^Sp and the adiabatic sound speed is Cg = P/p. Thus 
the total curvature perturbation is constant on large scales for purely adiabatic perturbations. 

In the presence of more than one fluid, the total non-adiabatic pressure perturbation, i5Pnad) niay be split into two 
parts. 


^Pnad — f^Pintr “b SPysI ■ 

The first part is due to the intrinsic entropy perturbation of each fluid 

'^Pintr — ^ ( ^Piiitr,a ; 


(27) 


(28) 


where the intrinsic non-adiabatic pressure perturbation of each fluid is given by 

<5Pintr.a = SPa - C^Spa , (29) 

c^ = Pal Pa is the adiabatic sound speed of that fluid and the total adiabatic sound speed is the weighted sum of the 
adiabatic sound speeds of the individual fluids. 


c 


2 

s 



(30) 


The second part of the non-adiabatic pressure perturbation m is due to the relative entropy perturbation between 
different fluids, denoted by Sap in Eq. 

^Prel = Qjj . 'y ( PaPp (c„ — SaP ■ 

^ a,/3 


( 31 ) 
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The time-dependence of the intrinsic entropy perturbation, 6Pinti,a, of each fluid must be specified according to 
the detailed modelling of that fluid. For instance, if the fluid has a definite equation of state = Pa{Pa) then the 
intrinsic non-adiabatic pressure perturbation vanishes. 

The evolution of the relative entropy perturbation, Saji, follows from the time dependence of the individual curvature 
perturbations Ca and Equation 11611 for the evolution of the gauge-dependent density perturbations in the long- 
wavelength limit reduces to 


^Pa + ^H{Spa + SPa) — 3 {pa + Pa) '*/’ + Qa4' + ^Qa ■ 


(32) 


Re-writing this in terms of the gauge-invariant curvature perturbation C,a defined in gives an evolution equation 
for the curvature perturbation on uniform a-fluid density hypersurfaces, 


Ca = [^Pa - “ 

Pa. 


HQa 

SQa 

( P _ 

Qa'^ 

5pa 

Pa 

Qa 

[2p 

Qa j 

Pa 


Po 


(33) 


For non-interacting, perfect fluids {Qa = 0 and i5Pintr,a = 0) we have Ca = 0 and the individual curvature perturbations 
for each fluid remain constant in the long-wavelength limit 0. But in general, the curvature perturbation ()„ may 
change with time either due to the intrinsic non-adiabatic pressure perturbation, dPintr.a in Eq. or due to what 
we will call the ‘non-adiabatic’ energy transfer, (5Qnad,a- 

Analogously to the total non-adiabatic pressure perturbation we will split the non-adiabatic energy transfer 
into two parts. 


^i5nad,a — T ^Qre\,a ■ 


(34) 


The first part is the instrinsic non-adiabatic energy transfer perturbations, defined as 


^Qintr.a = SQa “ ^5pa ■ 
Pa 


(35) 


This is automatically zero if the local energy transfer Qa is a function of the local density pa so that 5Qa = Qa^Paj Pa-, 
just as the intrinsic non-adiabatic pressure perturbation (HI vanishes when SPa = PaSpajpa- The second part is the 
relative non-adiabatic energy transfer 




QaP f ^Pa ^P 


Qa 

6 ^ 




a(3 


(36) 


where we have used the background Einstein equations 0 and the perturbed Friedmann constraint equation on 
large scales, 


^ + = (37) 

2 p 

in order to write (jQrei.a explicitly in terms of the relative entropy perturbation, Sap- 

Note that the relative non-adiabatic pressure perturbation, defined in Eq. m is related to the relative non-adiabatic 
energy transfer perturbation dehned in Eq. (TTHll as 

SPrel = -2^ V ^5Qrel,a • (38) 

P^ Qa 

The non-adiabatic pressure perturbations, (5Pintr,a and JPreh and the non-adiabatic energy transfers, ^Qintr.a and 
<5(3rei,a, are all automatically gauge-invariant. The intrinsic entropy perturbations, (5Pintr,a and (5Qintr,a, are both 
zero if the pressure and local energy transfer are determined by the local energy density. But even if the intrinsic 
entropy perturbations vanish, there may be a non-adiabatic energy transfer due to the relative entropy perturbation, 
C ~ Ca • We interpret this as due to a gravitational redshift (time-dilation) which perturbs the rate of energy transfer 
with respect to coordinate time if the uniform a-density hypersurface does not coincide with the uniform total density 
hypersurface, Ca ^ C- 
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By taking the difference between the evolution equations for two fluids we obtain an evolution equation for 
the relative entropy perturbation on large scales 


Sa/3 = 3iJ 


I 3IISPiiitT,a ^QintT,a J-Pintr,/? *^Qintr,/3 


P/S 


E P7 / o ^/s Q 

On n.f 


2p 


Qa Q/3 


PP 


Pi 


(39) 


Thus we see that any relative entropy perturbation Sap is sourced on large scales only by intrinsic entropy per¬ 
turbations in the a and (3 fluid, or by other relative entropy perturbations. There is no source term coming from 
the overall curvature perturbation and so adiabatic perturbations (with no intrinsic or relative entropy perturbation) 
remain adiabatic on large-scales even when one considers interacting fluids. 


III. CURVATON DECAY 

Having established a general formalism in which to study the evolution of large-scale curvature and entropy per¬ 
turbations including entropy transfer between multiple fluids, we now study the specific case of a non-relativistic 
matter decaying into radiation. In particular this can be used to describe the decay of a massive curvaton field into 
radiation [T^lift 1^1^. 

The curvaton scenario has recently been proposed 0, |3 1^ as a mechanism by which a large-scale curvature 

perturbation can be produced from an initially isocurvature perturbation. If the curvaton is a light scalar field (with 
mass less than the Hubble rate) the field may acquire an almost scale-invariant spectrum of perturbations, Co-. In 
the curvaton scenario radiation, p^, is supposed to dominate the initial energy density after inflation and this is 
assumed to be unperturbed, C 7 = 0. Thus the curvaton perturbation is initially an isocurvature density perturbation 
(C — 0 and Scr-y = 3Co) and remains an isocurvature perturbation while the relative density of the curvaton remains 
negligible. However, once the Hubble rate drops below the mass of the curvaton, the field begins to oscillate. Averaged 
over several oscillations the effective equation of state is {Pa/ Pa) = 0 , i.e., the coherent oscillations of the field are 
equivalent to a fluid of non-relativistic particles . As the energy density of non-relativistic particles grows relative 
to the energy density of radiation, what was once an isocurvature perturbation becomes a perturbation in the total 
curvature, Eq. ED- 

Assuming the curvaton is unstable and decays into light particles (“radiation”) with a decay rate T, this represents 
an energy transfer from the pressureless curvaton fluid to the radiation fluid. The precise amplitude of the resulting 
curvature perturbation, relative to the initial curvaton perturbation, depends upon both the initial density and 
the decay rate of the curvaton. We present the equations for the evolution of the curvature and relative entropy 
perturbations and solve them numerically, comparing with an analytic approximation assuming an instantaneous 
decay. 


A. Background solution 

The energy transfer from the massive curvaton to light radiation is described by 



Q<t — r P(j , 

(40) 


— r Pfj , 

(41) 

where T is the decay rate of the curvaton into radiation, which we take to be a constant. 

The energy conservation 

equations are therefore 

Pa = -Pa (3i? -1- r) , 

(42) 


p^ = -AHp^ + Tpa, 

(43) 


p = -H {3pa + ip^) , 

(44) 

where the Hubble expansion is given by 

2 SttG 

I? — 2 {Pa + P 7 ) • 

(45) 


In order to solve the system of equations above numerically, it is convenient to work in terms of the dimensionless 
density parameters 




( 46 ) 
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FIG. 1: 


Phase-plane showing trajectories for the background solutions in the curvaton model in Eqs. 1I48II51II . 


and the dimensionless “reduced” decay rate 

^ _ r 

^ ^ T + H ’ 

which varies monotonically from 0 to 1 in an expanding universe. 

The background equations 142I45II can then be written as an autonomous system 


(47) 


n' 


g' 



g)g, 


(48) 

(49) 

(50) 


where ' denotes differentiation with respect to the number of e-foldings N = lna. The density parameters are subject 
to the Friedmann constraint which requires 


= (51) 

There are only two independent dynamical equations and the generic solutions follow trajectories in a compact two- 
dimensional phase-plane (0 < 5 < 1, 0 < < 1), illustrated in Figure^ 

The dynamical system (I48H51II admits three fixed points 

(A) n^ = i,n, = o,g = 0 , 

(B) n^ = o^n^ = i,g = 0 , 

(C) = = g = 1 . 

Generic solutions start at the unstable repellor {A) and approach the stable attractor {C) at late times. At early 
times (ri-), ~ 1, 5 <C 1) we find 5 oc cx a^. The standard radiation dominated cosmology corresponds to evolution 
along the line rio- = 0. However solutions can approach arbitrarily close the curvaton-dominated saddle point {B) 
before the curvaton decays and Her —> 0 once again. 


B. Perturbations 

Both the curvaton and radiation fluids have fixed equations of state = 0 and 6P^ = Sp^/S) and hence there 
can be no intrinsic non-adiabatic pressure perturbation {SPintr.a = 0 and ( 5 Pintr ,7 = 0). However the total curvature 
perturbation, does change on large scales in the presence of a relative entropy perturbation (P|) 


5,^7 = 3(C. - C 7 ) : 


(52) 














which leads to a non-adiabatic pressure perturbation EU- The evolution of the total curvature perturbation using 
Eqs. (EHll . is 


c 


H Po'P'y 

y 


5. 


(53) 


We assume that the curvaton decay rate T is fixed by microphysics (i.e., JT = 0) and hence the perturbed energy 
transfer is simply given by 


6Qa = -rSpa, (54) 

5Qj = TSpa- ■ (55) 

This energy transfer is determined solely by the local density of the curvaton and hence there is no intrinsic non- 
adiabatic energy transfer from the curvaton, i^Qintr.cr = 0. However the radiation suffers an intrinsically non-adiabatic 
energy transfer from the curvaton decay ll^ 


'5Qintr,7 = T [5pa- ^^Pl \ > 


(56) 


which is proportional to the relative entropy perturbation between the radiation and curvaton 


^Qintr,7 — cy jt ’ 

orl 


(57) 


The relative non-adiabatic energy transfers EHl) are also non-zero and given by, 


^Crel,cr 


(5Qrel, 


7 


^PaP ( Sp^ _ 

‘^P \ Pa P ) 

^PaP f ^ 

2p V P / ’ 


which can be rewritten in terms of the relative entropy perturbation as 


*^Crel,{7 

f^Crel,7 


^Pa 

6Hp 

6Hp 


PjSa^ , 

Pa^a-y ■ 


(58) 

(59) 


(60) 

(61) 


Thus the evolution equations for the curvature perturbation on uniform curvaton density hypersurfaces, C,a-, and 
uniform radiation density hypersurfaces, are given by 


C 7 


_r pct P 7 

„ . <Jay ) 

& P Pa 

r Pa C _ Pa \ 

3 P 7 V 2p) 


Sa 


(62) 

(63) 


The evolution equation for the relative entropy perturbation Say is, from Eq. (I39II . 


b Pa Pa 

2 Py P 



Sa 


(64) 


Equations 15311 and (IMIl form a closed system of first-order equations for the evolution of the adiabatic and entropy 
perturbations, and Sap, on large-scales in the curvaton model. They clearly demonstrate the general principle 
that the total curvature perturbation, Q, evolves on large scales in the presence of a relative entropy perturbation, 
Sap, while the entropy perturbation obeys a homogeneous evolution equation, unaffected by the large-scale curvature 
perturbation. Alternatively we could use Eqs. (EH) and Ol as a closed system of first-order equations for (^a and C 7 , 
remembering that Say = ^{Ca — Cy)- 

However the evolution equation ESJ for and, hence, the evolution equation EH) for Say both become singular 
whenever Tpa = iHp^ and p^ = Q. This is due to the uniform p^ hypersurface becoming ill-defined rather than any 
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breakdown of perturbation theory on generic hypersurfaces. In particular the uniform-po- and uniform total energy 
density hypersurfaces remain well-behaved. 

In practice we will use the two non-singular evolution equations (I5d|l and (PI for C and (a-, respectively. In terms 
of the dimensionless background variables we have two coupled evolution equations 


C' 


c 


Oa(2g-3) 
g(4 - 


(C - Ca), 


2(3 - 2g) 


(C - C^) • 


(65) 

( 66 ) 


To calculate the final curvature perturbation on large scales produced in the curvaton scenario we start with initial 
conditions close to the point (A) for the background variables (g <C 1 , Oo- <C 1 ) and unperturbed radiation, but 
perturbed curvaton fluid: 


C 7 = 0, (67) 

G = Ca.in- (68) 

From the definitions of the total curvature perturbation and the entropy perturbation Eqs. and P|). this corre¬ 
sponds to initial values 


c 


■^(77 


■ Ccr.in ) 


3n. 


4-0. 

3C^,in. 


(69) 

(70) 


This is an initial isocurvature perturbation in the sense that C 0 in the early time limit, (A) , where > 0. 

Starting from these initial conditions we use Eqs. (PI and to follow the evolution of C and until we reach 
the late time attractor (C) where g ^ I and O. ^ 0. At late times the perturbations too approach a fixed point 
attractor where and 


C — ^Co-.in , (71) 

5,.^ = 0. (72) 

This is an adiabatic primordial perturbation, where the final value of the large-scale curvature perturbation, C, is 
related to the initital curvaton perturbation Co-,in by a parameter r MM which is determined by the numerical 
solution of Eqs. (ESI) . itSni) . (IS51) and (ESI) . 

Thus, we can represent the integrated effect upon the large-scale curvature and entropy perturbations of the curvaton 
growth and decay by the transfer matrix: 


c 


5, 


(77 



1 r/3 

0 0 



(73) 


Examples of the evolution of large-scale perturbations for two different choices of initial conditions are shown in 
Figures |21 and |21 The resulting value for the transfer parameter r defined in Eq. GD depends upon the maximum 
value of fla- before the curvaton decays. If the curvaton dominates before it decays, i.e., 12.,dec — 1, we have r ~ 1 as 
in the case shown in FigEl More generally, r is a one-dimensional function of the initial value of f2./(F/i7)^/^ which 
determines which trajectory is followed in the two-dimensional (g, 12.) phase-plane, Figure^ The precise dependence 
of r upon the initial value of f2./(F/i7)^/^ is shown in Figure^] 


C. Comparison with sudden decay approximation 

Previous analyses [Sllllli have relied on the assumpton of “sudden decay” to estimate the final curvature 
perturbation produced after curvaton decay. In this approximation the energy transfer (5. = —Fp. is assumed to be 
negligible until the decay time, defined by F /H reaching some critical value F /H^ec of order unity, at which time all 
the energy density of the curvaton field is rapidly converted into radiation. 

In the absence of energy transfer the individual curvature perturbations C. and C 7 ) defined by Eq. (EUI) , remain 
constant on large scales [see Eq. (I 2 SI]- Thus the total curvature perturbation, Eq. is given by 


C « /C. + (1 - /)C7, 


(74) 
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FIG. 2: Evolution of the normalised curvature perturbation on uniform curvaton density hypersurfaces, ^c/Co-.in, and of the 
normalised total curvature perturbation, C/(^CT,in, as a function of the number of e-foldingss, starting with Co-/C<T,in = 1 and 
initial density and decay rate = 10“^ and T/H = 10“®. 



FIG. 3: Evolution of the normalised curvature perturbation on uniform curvaton density hypersurfaces, ^cr/Co-.in, and of the 
normalised total curvature perturbation, ^/(^<T,in, as a function of the number of e-foldingss, starting with ^o-/C<T,in = 1 and 
initial density and decay rate 0,^ = 10“^ and T/H = 10“®. 



FIG. 4: Transfer parameter r defined in Eq. ffn obtained from numerical solutions as a function of the initial value of 
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where Co- and C7 are constant and the only time dependence arises from the time dependence of the weight given to 
the curvaton perturbation 


After the curvaton decays into radiation all the energy density in the model has a unique equation of state, hence 
SPna.d = 0 in Eq. 112611 ■ and ( becomes constant on large scales. Hence in the curvaton scenario, where the initial 
curvature perturbation is assumed to be negligible, the resulting adiabatic curvature perturbation after curvaton 
decay is given by 


Cout ~ fdecC(7,in • 

In terms of an initial value for Ho-,in and reduced decay rate r/iJjn, we can write /dec as 


/dec — 


3H„,in 


“1” 4(1 Hn'dn)ydec 


(76) 


(77) 


where ydec = «in/adec is the ratio of initial scale factor to that at decay. We can calculate this from the Friedmann 
equation for non-interacting matter and radiation which can be written as 






Thus the epoch of decay, j/deo is given by the one real root, 0 < j/dec < 1 of 


(1 ^o-,in)ydec 4“ i^n-jn Udec 


( r/i7ir 


Vr/i7de 


= 0 , 


and /dec is then obtained from Eq. iTTTI) . 
There are two limiting cases: 


J - f^a,in) ( y/gr ) ^dec > ^^<T,in/(l “ f^fT.in) 

/dec ~ \ , / , \ 1/2 

1 - (1 - ^^-.in) for 2/dec « f^..in/(l " f^^.in) 


(78) 


(79) 


(80) 


The sudden decay approximation is compared with numerical results for the full equations (ESll and in Figure 0 
The one free parameter in the sudden-decay approximation is the particular value of T/TJdec chosen to characterise 
the epoch of decay. Optimising the fit to the full numerical solutions fixes T/Hdec — 1-4, in line with our expectation 
that r/i7dec should be of order unity. With this choice the sudden-decay approximation is seen (Figure [SJ to give a 
good estimate for r ~ /dec (good to within 10%). 


IV. CONCLUSIONS 

We have studied the evolution of large-scale curvature perturbations for multiple interacting fluids in a linearly per¬ 
turbed FRW cosmology. The curvature perturbation, on hypersurfaces of uniform density for each fluid, Eq. CT . 
provides a gauge-invariant variable by which to study the large-scale evolution. The total curvature perturbation, 
is then a weighted sum, Eq. of the individual /a’s. For a non-interacting perfect fluid, C,a remains constant on 
large scales, independently of perturbations in other fluids 0- More generally we have shown how C,a can change on 
large scales due to either an intrinsic non-adiabatic pressure perturbation or non-adiabatic energy transfer. 

We can decompose an arbitrary energy transfer perturbation into two parts: 

5Qa = ^Spa + SQa, intr , (81) 

Pol 

where 5Qa,intr is the (gauge-invariant) intrinsic non-adiabatic energy transfer. The large-scale curvature perturbation 
Qa can change due to this intrinsic non-adiabatic energy transfer, or due to a relative entropy perturbation between 
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FIG. 5: Comparison of full numerical solution for r in Eq. IZHJ with sudden-decay approximation, /dec given in Eq. IZZJ. 


fluids, (JQa.rei defined in Eq. which is proportional to Ca ~ C- For perturbations that obey the generalised 

adiabatic condition: 

SPa = c\5pa , 5Qa = and Ca = C, (82) 

Pa 

the curvature perturbation, remains constant on large scales. If all the individual fluids obey this generalised 
adiabatic condition, then the total curvature perturbation, C, is necessarily constant too. 

Most previous analyses have adopted the variables of Kodama and Sasaki who defined the relative entropy 
perturbation between two fluids as 


_ 5 Pa 5pp 

I D I D ' 

Pa H" J^a Pp + 


(83) 


However this definition is only gauge-invariant in the absence of energy transfer. As a result the evolution equations 
including energy transfer are particularly unpleasant 1^ . In particular it is diffcult to show that entropy pertur¬ 
bations obey a homogeneous evolution equation on large scales (i.e., adiabatic pert urbations stay adiabatic on large 
scales) in the way that was recently shown for multiple interacting scalar fields Although it has been be argued 
on very general grounds that this must be the case 0, this fundamental result has not previously been explicitly 
demonstrated in treatments of interacting fluids. 

We have used the correct gauge-invariant generalisation of allowing for energy transfer. 


Sap = 3(Ca - C/?): 


(84) 


which describes the relative displacement between the two hypersurfaces of uniform density defined with respect to 
the two fluids. This reduces to 1I83|1 in the case of no energy transfer. It allows us to demonstrate that the evolution of 
the large-scale entropy perturbation, Eq. is sourced only by entropy perturbations and not sourced by the total 
curvature perturbation, Q. Thus the integrated evolution on large scales, even when we include energy transfer, can 
be schematically represented by the linear transfer matrix: 


c 

5 



1 T^s 
0 Tss 



(85) 


We have applied our formalism to study the evolution of curvature perturbations in the curvaton scenario where 
an initially isocurvature (non-adiabatic) perturbation in the curvaton field is transferred to the radiation fluid when 
the curvaton eventually decays. The decay of the curvaton represents a non-adiabatic energy source for the radiation 
fluid. We have numerically solved the coupled evolution equations to determine the resulting curvature perturbation, 
C, for an initial entropy perturbation. Thus we have calculated the transfer coefficient Tqs in Eq. 18511 for different 
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parameter values of the background models. We compared our results with semi-analytic estimates based on the 
“sudden-decay” approximation [l^l^ where the fluids are assumed to be non-interacting up until a fixed decay time. 
The sudden-decay approximation is shown to give a good fit to the full result (within 10%) for a suitable choice of 
fitting parameter. 

In this two-fluid realisation of the curvaton scenario the interaction between the fluids leads to the relative entropy 
decaying to zero at late times, Tss = 0 in Eq. leaving a purely adiabatic curvature perturbation. Our formalism 
can also be applied to cosmological models including other cosmological fluids such as baryons, CDM or neutrinos, 
in which case it should be possible to calculate the amplitude of any residual isocurvature perturbations that may 
survive after curvaton decay in different variations of the curvaton scenario. 
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